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The mixed spin-1/2 and spin-3/2 Ising model on the extended Kagome lat- 
tice is solved by establishing a mapping correspondence with the eight-vertex 
model. When the parameter of uniaxial single-ion anisotropy tends to infinity, 
the model system becomes exactly soluble as the staggered eight-vertex model 
satisfying the free-fermion condition. The critical points within this manifold 
can be characterized by critical exponents from the standard Ising universality 
class. The critical points within another subspace of interaction parameters, 
which corresponds to a coexistence surface between two ordered phases, can 
be approximated by corresponding results of the uniform eight-vertex model 
satisfying the zero-field condition. This coexistence surface is bounded by a line 
of bicritical points that have non-universal continuously varying critical indices. 
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1. Introduction 

Investigation of phase transitions and critical phenomena belongs to the most in- 
tensively studied topics in the equilibrium statistical physics. A considerable progress 
in the understanding of order- disorder phenomena has been achieved by solving 
planar Ising models that represent valuable exceptions of exactly soluble lattice- 
statistical models with a non-trivial critical behaviour [ 1 . Although phase transi- 
tions of planar Ising models have already been understood in many respects, there 
are still a lot of obscurities connected with a criticality of more complicated spin 
systems exhibiting reentrant transitions, non-universal critical behaviour, tricritical 
phenomenon, etc. It is worthy to mention, however, that several complicated Ising 
models can exactly be treated by transforming them to the solvable vertex models. A 
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spin-1/2 Ising model on the union jack (centered square) lattice, which represents a 
first exactly soluble system exhibiting reentrant transitions ['SI, can be for instance 
reformulated as a free-fermion eight- vertex model [jSj. It should be also pointed out 
that an equivalence with the vertex models have already provided a precise con- 
firmation of the reentrant phenomenon in the anisotropic spin-1/2 Ising models on 
extended Kagome lattice [ E] and centered honeycomb lattice [ 5J as well. 

Despite the significant amount of effort, there are only few exactly soluble Ising 
models consisting of mixed spins of different magnitudes, which are usually called 
also as mixed-spin Ising models. A strong scientific interest focused on the mixed- 
spin systems arises partly on account of much richer critical behaviour they display 
compared with their single-spin counterparts and partly due to the fact that they 
represent the most simple models of ferrimagnets having a wide potential applicabil- 
ity in practice. Using the extended versions of decoration-iteration and star-triangle 
mapping transformations. Fisher [IHj has derived exact solutions of the mixed spin- 
1/2 and spin-S* {S > 1) Ising models on the honeycomb, diced and decorated hon- 
eycomb lattices. Notice that these mapping transformations were later on further 
generalized in order to account also for the single-ion anisotropy effect. The influence 
of uniaxial and biaxial single-ion anisotropies have precisely been investigated on the 
mixed-spin honeycomb lattice [13 as well as on some decorated planar lattices [IH]. 
With exception of several mixed-spin models formulated on the Bethe (Cayley tree) 
lattices, which can be accurately treated within a discrete non-linear map [ IT or 
an approach based on exact recursion equations [^^, these are the only mixed-spin 
planar Ising models with generally known exact solutions, yet. 

One of the most outstanding findings emerging in the phase transition theory 
is being a non-universal critical behaviour of some planar Ising models, which is in 
obvious contradiction with the idea of universality hypothesis [E]. The mixed spin- 
1/2 and spin-S* Ising model on the union jack lattice [E] represents very interesting 
system from this viewpoint as it exhibits a remarkable line of bicritical points that 
have continuously varying critical indices obeying the weak universality hypothesis [ 
IT^ . In the present article, we shall investigate a topologically similar mixed spin-1/2 
and spin-3/2 Ising model on the extended Kagome lattice by establishing a mapping 
correspondence with the staggered and uniform eight-vertex models, respectively. In 
a certain subspace of interaction parameters, the model under investigation becomes 
exactly soluble as the staggered eight- vertex model satisfying the free-fermion condi- 
tion [|14j. Even if a non- validity of the free-fermion condition in the rest of parameter 
space is simply ignored, one still obtains rather reliable estimate of the criticality 
within free-fermion approximation [E]- Finally, the critical points within another 
subspace of interaction parameters can be approximated from the relevant solution 
of the uniform eight- vertex model satisfying the zero-field condition. 

The outline of this paper is as follows. In Section 2, a detailed formulation of 
the model is presented and subsequently, the mapping correspondence that ensures 
an equivalence with the eight-vertex models will be derived. The most interesting 
numerical results for a critical behaviour will be presented and particularly discussed 
in Section 3. Finally, some concluding remarks are drawn in Section 4. 
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Figure 1. Diagrammatic representation of the extended Kagome lattice com- 
posed of the mixed spin-1/2 (empty circles) and spin-3/2 (filled circles) sites, 
respectively. The solid (broken) lines depict the nearest-neighbour (next-nearest- 
neighbour) interactions. 



2. Formulation 

Let us begin by considering the mixed spin-1/2 and spin-3/2 Ising model on 
the extended Kagome lattice C schematically illustrated in figure H The mixed-spin 
Kagome lattice consists of the spin-1/2 (empty) and spin-3/2 (filled circles) atoms 
placed on the six- and four-coordinated sites, respectively. The total Hamiltonian 
defined upon the underlying lattice C reads: 

2N 2N N/2 

'Hmix = -J X! SiCTj-J' ^ akai-D^Sl, (1) 

ik,l)ClC 1=1 

where aj = ±1/2 and Si = ±1/2, ±3/2 are Ising spin variables, J denotes the 
exchange interaction between nearest-neighbouring spin-1/2 and spin-3/2 pairs and 
J' labels the interaction between the spin-1/2 pairs that are next-nearest-neighbours 
on the extended Kagome lattice C. Finally, the parameter D measures a strength 
of the uniaxial single- ion anisotropy acting on the spin-3/2 sites and N denotes the 
total number of the spin-1/2 sites. 

In order to proceed further with calculation, the central spin-3/2 atoms should be 
firstly decimated from all faces of extended Kagome lattice After the decimation, 
i.e. after performing a summation over spin degrees of freedom of the spin-3/2 sites 
(filled circles), the partition function of the mixed-spin system can be rewritten as: 

N/2 N/2 

Zmix = XI n ^m(^j5 0-j,o-fc,cr/) J| uj^{ai, a j , ak, ai) . (2) 

{tr} m=l n=l 

Above, the summation is performed over all possible spin configurations available at 
the spin-1/2 sites and the first (second) product is over A^/2 faces having four spin- 
1/2 sites (jj, (jj, (jfc, ai placed in the corners of square plaquettes with (without) a 
central spin-3/2 site in the middle of these plaquettes (see figure H]). The Boltzmann 
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Figure 2. The eight possible line arrangements at a vertex of the dual lattice. 

factors to'^{a, b, c, d) and uj'^{a, b, c, d) assigned to two different kinds of alternating 
faces, which constitute the checkerboard lattice, can be defined as: 

uj-^{a,b,c,d) = 2exp[K'{ab + bc + cd + da)/2 + A/A] 

{exp(2A) cosh[3if (a + b + c + d)/2] + cosh[ir(a + b + c + d)/2]}, 
uj^{a,b,c,d) = exp[K'{ab + be + cd + da)/2], (3) 

where K = J/{k-QT), K' = J'/lk^T), A = D/{kBT), kB is Boltzmann's constant, 
and T stands for the absolute temperature. 

At this stage, the model under investigation can be rather straightforwardly 
mapped onto the staggered eight-vertex model defined on a dual checkerboard lat- 
tice Ct), since Boltzmann factors uj-^{a,b,c,d) and u^{a,b,c,d) are being invariant 
under the reversal of all four spin variables. Actually, there are maximally eight 
different spin arrangements giving different Boltzmann weights iv^{a, b, c, d) and 
u;^(a, b, c, d) for each kind of face. Diagrammatic representation of eight possible 
spin arrangements and their corresponding line coverings of the eight- vertex model 
is shown in figure El If, and only if, the adjacent spins are aligned opposite to each 
other, then solid lines are drawn on the edges of the dual lattice £7?, otherwise 
they are drawn as broken lines. It can be easily understood that eight possible line 
coverings around each vertex of the dual checkerboard lattice always correspond to 
two spin configurations, one is being obtained from the other by reversing all spins. 
Since there is even number of solid (broken) lines incident to each vertex of the dual 
lattice Cv, the model becomes equivalent to the staggered eight-vertex model. 

The Boltzmann weights uj-^{a, b, c, d) and Lj^{a, b, c, d), which correspond to eight 
possible line coverings emerging at vertices of the dual checkerboard lattice, can 
directly be calculated from equation Q: 

uj^ = 2exp(K72 + A/4)[exp(2A) cosh(3X) + cosh(ir)], 
cu^ = 2exp{~K'/2 + A/4)[exp(2A) + 1], 

= cu^ = 2exp(A/4)[exp(2A) + 1], 
ujf = uj^ = uj:^ = uj^ = 2exp(A/4)[exp(2A) cosh(3if/2) + cosh(ir/2)]; (4) 
ujf = exp{K'/2), cjf = exp{~K'/2), 

ujf = cuf = cuf = UJ^ = uo^ = Loi = 1. (5) 

Unfortunately, there does not exist general exact solution of the staggered eight- 
vertex model with arbitrary Boltzmann weights cof' and to^ = 1—8). However, 
if the weights ^ and (jSJ satisfy so-called free-fermion condition: 

^1^2 + ^Is^i = ^5^l6 + ^7^8, (6) 
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the staggered eight-vertex model then becomes exactly soluble as the free-fermion 
model solved several years ago by Hsue, Lin and Wu [UM. The expressions which 
enter into the free-fermion condition © can be defined through: 

00 - , A A -u , A, ,A ,^ -u , A fi, A, P -u , A, P, A 

iiSiig — ^\ ^3 ^3 I ^2 ^2 ^4 "^4 ~r ^5 ^6 ^7 ^8 ^5 ^6 7 8 ' 

It can be readily proved that the free-fermion condition © holds in our case just as 
D— >-±oo, orT— i>oo. The restriction to infinitely strong single-ion anisotropy con- 
sequently leads to the familiar phase transitions from the standard Ising universality 
class, because in this case our model effectively reduces to a simple spin-1/2 Ising 
model on the extended Kagome lattice. Within the manifold given by the constraint 
(jH)), the free-fermion model becomes critical as long as: 

VLi + VL2 + VL3 + VLi = 2max{fii, VL2, ^3, ^a}- (8) 

It is noteworthy, however, that the critical condition (jSJ yields rather reliable es- 
timate of the criticality within so-called free-fermion approximation [ even if a 
non-validity of the free-fermion condition © is simply ignored. 

Now, we shall establish an approximate mapping between the staggered and 
uniform eight-vertex models, since the second branch of exact solution is available 
just for the latter model under the zero-field condition [[T]. For this purpose, let us 
define average Boltzmann weights of the staggered eight- vertex model, which would 
approximately transform the staggered eight-vertex model into the uniform one: 

u, = u:fu:f, (z = l-8). (9) 

Note that the uniform eight- vertex model satisfies the zero-field condition just when 
its Boltzmann weights are pairwise and symmetrically equal to each other: 

Ui = ^2, UJ3 = t^4, = ^6, ^7 = ^8- (10) 

As we already have = 0)4, 0)5 = uq, and 0)7 = cDg, the zero-field case is consequently 
reached by imposing the condition uji = 0)2 only, or equivalently: 

, , , exp(-2K') - coshffs:) , , 

cosh(3A ) — exp(— 2a') 

According to Baxter's exact solution [ilj, the zero-field eight- vertex model becomes 
critical on the manifold pO|) if: 

uji + UJ3 + Ojr^ + uj-j = 2max{cJi, Oj^}. (12) 

It is easy to check that uji represents in our case the largest Boltzmann weight, thus, 
the condition determining the criticality can also be written in this equivalent form: 

exp{K'^) [exp(2AJ cosh(3/s:j + cosh{K^)] = 

1 + exp(2Ac) + 2exp(2Ac) cosh(3is:e/2) + 2 cosh(irc/2), (13) 
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Figure 3. Ground-state phase diagram in the J' — D plane when J > 0. Broken 
rectangles schematically illustrate a typical spin configuration within each phase. 
Broken line connecting both triple points shows a projection of the critical line 
((131) into the J' - D plane. 
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where K^. = J/^k-QTc), K'^ = J' /{k-oTc), = D/ik-oTc), and denotes the critical 
temperature. It should be stressed, nevertheless, that the critical exponents (with 
exception of 5 and rj) describing a phase transition of the zero-field eight-vertex 
model depend on the function ji = 2 arctan(a;5a;7/(I;ia;3)^/^, in fact: 

a = a' = 2--, P=T^, u = u' = ^, 7 = ^, ^ = 15, V = ], (14) 

Finally, let us explicitly evaluate the critical exponent P that determines disappear- 
ance of the spontaneous order as the critical temperature is approached from below: 

1^32 f exp(2A,) cosh(3ire/2) + cosh(ir,/2) \ 

IJ vr ^^""^^^t [exp(2A,) + l]3/4[exp(2A,) cosh(3i^e) + cosh(ir,)] V4 f ^ ^> 



3. Results and discussion 



Now, let us turn our attention to a discussion of the most interesting results 
obtained for the ground-state and finite-temperature phase diagrams. Solid lines 
displayed in figureElrepresent ground-state phase boundaries separating four distinct 
long-range ordered phases that emerge in the ground state when J > 0. Spin order 
drawn in broken rectangles shows a typical spin configuration within basic unit cell 
of each phase. As could be expected, a sufficiently strong antiferromagnetic next- 
nearest-neighbour interaction J' alters the structure of the ground state owing to 
a competing effect with the nearest-neighbour interaction J. Due to a competition 
between the interactions, the central spins are free to flip within the phases III and 
IV and thus, these phases exhibit a remarkable coexistence of order and disorder. 
At last, it is worthwhile to mention that a broken line connecting both triple points 
depicts a projection of the approximate critical line ()13|1 into the J' — D plane. As 
this projection crosses zero-temperature plane along the ground-state transition line 
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Figure 4. (a) The curve scaled to the left axis shows how the critical temperature 
changes with a strength of the next-nearest-neighbour coupling J'/ J, the curve 
scaled to the right axis depicts a variation the single-ion anisotropy parameter 
along this line; (b) The same as for figureEta), but the critical exponent (3 is now 
scaled to the left axis. Broken lines are in both figures guides for eyes only. 



D/ J = —3/2 — J' / J between the phases I and IV, it is quite reasonable to suspect 
that this line represent a location of phase transitions between these phases. 

Let us investigate more deeply this line of critical points. The critical tempera- 
tures calculated from the uniform zero-field eight- vertex model must simultaneously 
obey both the zero-field condition (|TT|l as well as the critical condition ()13j) . It 
is easy to check that the former condition necessitates —1.5 < J' /J < —0.5 and 
— 1.0 < D/J < 0.0. Figure IHa) displays a projection of this critical line into the 
J' — Tc plane (the dependence scaled to the left axis) and respectively, a projection 
into the J' — D plane which is scaled to the right axis. Along this critical line, the 
critical exponents are expected to vary with interaction parameters as they have to 
follow the equations For illustration, figure Efb) shows how the critical index (3 
changes along the critical line. Apparently, the exponent (3 approaches its smallest 
possible value 1/16 by reaching both triple points with zero critical temperature, 
however, it is also quite interesting to ascertain that its greatest value is below the 
value 1/8 that predicts the universality hypothesis for planar Ising systems [C]. 

Before concluding, few remarks should be addressed to a global finite-temperature 
phase diagram plotted in figure El which displays the critical temperature as a func- 
tion of the ratio J' / J for several values of the single-ion anisotropy D/J. Critical 
boundaries depicted as solid lines represent exact critical points obtained from the 
free-fermion solution (jH)) of the staggered eight- vertex model obtained under the con- 
straint which is fulfilled in the limiting cases D/J — ±00. Dotted critical lines 
show estimated critical temperatures calculated from the free-fermion approximation 
simply ignoring a non-validity of the free-fermion condition © for any finite value 
of D/J. Approximative solution related to the critical points (fT^ of the uniform 
zero-field eight-vertex model on the variety (fTT|) is displayed as a rounded broken 
line. It is quite obvious from the ground-state phase diagram (figure that a right 
(left) wing of the displayed critical boundaries corresponds to the phase I (III) if 
D/J > 0.0, while it corresponds to the phase II (IV) if D/J < —1.0. Actually, the 
exact as well as approximate critical points resulting from the free-fermion solution 
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Figure 5. A plot critical temperature versus the ratio J'/ J for several values of 
the single-ion anisotropy term Dj J . For details see the text. 



correctly reproduce the ground-state boundaries between these phases. When the 
single-ion anisotropy parameter is selected within the range —1.0 < Dj J < 0.0 (see 
for instance the curve for D/J = —0.5), however, the critical line obtained from 
the free-fermion approximation meets at a bicritical (circled) point with the critical 
line of the equivalent uniform zero-field eight-vertex model as it has been already 
reasoned by Lipowski and Horiguchi [^TT who have solved similar spin system on 
the union-jack lattice. In such a case, the right and left part (with respect to the 
bicritical point) of this critical line separate the phases I and IV, respectively, and 
a line of first-order phase transitions is expected to terminate at this special mul- 
ticritical point. There are strong indications supporting this concept [lEI, actually, 
the almost straight broken line depicting the zero-field condition should always 
show a coexistence of these two phases as it starts from a point that determines their 
coexistence in the ground state. With regard to the aforementioned arguments one 
may conclude that a coexistence surface between the phases I and IV lies inside the 
area, which is bounded by the line of bicritical points (rounded broken line) having 
the non-universal interaction-dependent critical exponents. 

Finally, we should remark a feasible appearance of reentrant transitions which 
can be observed in the critical lines nearby the coexistence points D/J = 0.0 and 
—1.0. It is quite apparent that the observed reentrance can be explained in terms 
of the coexistence of a partial order and partial disorder emerging in both the high- 
temperature reentrant phases III and IV. As a matter of fact, the partial disorder 
of the spin-3/2 atoms can compensate a loss of entropy that occurs in these phases 
due to a thermally induced partial ordering of the spin-1/2 atoms what is in a good 
accordance with a necessary condition conjectured for the appearance of reentrant 
phase transitions 



4. Concluding Remarks 

The work reported in the present article provides a relatively precise information 
on the critical behaviour of the mixed spin-1/2 and spin-3/2 Ising model on the ex- 
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tended Kagome lattice by establishing a mapping correspondence with the staggered 
and uniform eight-vertex models, respectively. The main focus of the present work 
has been aimed at the examination of the criticality depending on the single-ion 
anisotropy strength as well as the strength of the competing next-nearest-neighbour 
interaction. The location of the critical boundaries has accurately been determined 
from the free-fermion solution of the staggered eight- vertex model and the zero-field 
solution of the uniform eight- vertex model, respectively, whereas the validity of both 
mappings is restricted to the certain subspaces of interaction parameters only. In the 
rest of parameter space, the free-fermion approximation has been used to estimate 
the critical boundaries as this method should provide meaningful approximation 
giving rather reliable estimate to the true transition temperatures. 

The greatest theoretical interest in this model arises due to the remarkable crit- 
ical line consisting of bicritical points, which bounds a coexistence surface between 
two long-range ordered phases. The bicritical points can be characterized by non- 
universal interaction-dependent critical exponents that satisfy the weak universal- 
ity hypothesis. Moreover, the same arguments as those suggested by Lipowski and 
Horiguchi have enabled us to identify the zero-field condition (jllj) with a loca- 
tion of the first-order transition lines separating these two ordered phases. 

It should be remarked that the considered spin system also shows reentrant 
phase transitions on account of the competition between the nearest- and next- 
nearest-neighbour interactions. Our results are in agreement with the conjecture [E] 
stating that the reentrance appears as a consequence of the coexistence of a partial 
order and disorder, namely, the partial disorder induced among spin-3/2 atoms can 
compensate the loss of entropy, which occurs on behalf of the partial ordering of the 
spin-1/2 atoms in both the high-temperature partially ordered phases. 
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